Single Cooper-pair pumping in the adiabatic limit and beyond 
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We demonstrate controlled pumping of Cooper pairs down to the level of a single pair per cycle, 
using an rf-driven Cooper-pair sluice. We also investigate the breakdown of the adiabatic dynamics 
in two different ways. By transferring many Cooper pairs at a time, we observe a crossover between 
pure Cooper-pair and mixed Cooper-pair-quasiparticle transport. By tuning the Josephson coupling 
that governs Cooper-pair tunneling, we characterize Landau-Zener transitions in our device. Our 
data are quantitatively accounted for by a simple model including decoherence effects. 



Charge pumps [IH3] and turnstiles [5] have recently 
attracted considerable attention. They could be used 
as building blocks for quantum computing devices [B], 
or to create a quantized current source that would pair 
up with the Josephson voltage and quantum Hall resis- 
tance to close the so-called quantum metrology triangle 
[7]. Among different types of realizations, Cooper-pair 
pumps [5HTU] stand out as macroscopically coherent ob- 
jects, with the phase of the superconducting order param- 
eter in the leads playing a key role. In addition, the cyclic 
path described in the space of parameters when pumping 
is equipped with a nontrivial geometric structure, allow- 
ing for the observation of geometric-phase effects [IT] . In 
the adiabatic limit, a general relation was derived [12 
connecting the pumped charge to the geometric (Berry) 
phase accumulated by the system ground state along a 
pumping cycle. This relation was experimentally demon- 
strated in Ref. 13, Beyond the adiabatic limit, we have 
recently proposed to employ a Cooper-pair pump as a 
Landau-Zener interferometer for geometric phases [*H] . 

In this Letter, we demonstrate controlled pumping of a 
single Cooper pair, using an rf-driven Cooper-pair sluice 
[HI ES] . Accessing this regime opens new possibilites for 
Cooper-pair pumping, from quantum metrology to the 
study of dissipation in driven quantum systems |16j , also 
in connection with geometric phases. We then investi- 
gate the breakdown of adiabatic pumping. In the sluice, 
this is expected to take place via Landau-Zener transi- 
tions (LZTs) at level anti-crossings. We reach the nona- 
diabatic limit in two different ways. By pumping many 
Cooper pairs at a time, we witness a crossover between 
pure Cooper-pair and mixed Cooper-pair-quasiparticle 
dynamics, due to continuous generation of nonequilib- 
rium quasiparticles by the nonadiabatic drive. By tun- 
ing the Josephson coupling that governs Cooper-pair tun- 
nelling, we characterize LZTs in our device. Our data are 
quantitatively accounted for by a simple model compris- 
ing LZTs and realistic decoherence. 

The Cooper-pair sluice (Fig. [TJ is a fully tunable 
Cooper-pair transistor, consisting of a small supercon- 
ducting island connected to leads by two superconducting 
quantum interference devices (SQUIDs). The SQUIDs 
are controlled independently by inductively coupled on- 




FIG. 1. (Color online) The Cooper-pair sluice, (a) Equivalent 
circuit of the sluice and scheme of the measurement set-up. 
(b) False-color scanning-electron micrograph of a representa- 
tive device, (c) Time evolution of the control parameters for 
a typical pumping cycle. 

chip coils, so that they can serve as Josephson junc- 
tions of tunable energy J%, J 2 . A gate electrode ca- 
pacitively coupled to the island induces a polarization 
charge n g = C g V g /2e in units of Cooper pairs, where V g 
is the gate voltage and C g the cross-capacitance between 
the gate and the island. The device is typically oper- 
ated in the charging regime, meaning that AEc S> J\, J 2 
(Ec — e 2 /2Cs is the single-electron charging energy of 
the island, Cs being its total capacitance). Pumping is 
realized by steering the three control parameters J\, J2 
and n g in a periodic fashion. The gate is used as a pis- 
ton to change the number of Cooper-pairs on the island, 
while the SQUIDs are operated as valves so as to impart 
a direction to the flow of charge. A typical pumping cycle 
is described in Fig.[ljc). 

We fabricate the devices by standard electron-beam 
lithography, two-angle Al evaporation, and liftoff. Small 
Josephson junctions (area ~ 70x70nm) are obtained by 
oxidization of the first Al layer in controlled O2 atmo- 
sphere. A scanning-electron micrograph of a representa- 
tive device is shown in Fig. [TJb). From the Coulomb- 
blockade conductance peak measured at 2 K we obtain 
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FIG. 2. (Color online) Single Cooper-pair pumping. Pumped 
charge Q p versus peak-to-peak amplitude An g and offset n° g 
of the gate drive. Dashed lines enclose the diamond-shaped 
regions where Q p /2e is expected to be constant and quantized 
(assuming ideal operation). 



Eq = 0.77 K. The normal-state resistance of the de- 
vice at 2K is Rn — 29kf2. Using the measured su- 
perconducting gap at base temperature A = 180 ueV, 
Rn and the Ambegaokar-Baratoff formula, we estimate a 
maximum Josephson energy J® lax — 0.46 K per SQUID. 
From switching statistics to the normal state in a current- 
biased configuration, we estimate the ratio between the 
maximum and minimum Josephson couplings obtained 
by varying the flux to be J m inl^max — 0-03 for both 
SQUIDs. 

All measurements are performed in a dilution refrigera- 
tor down to 20 mK. The set-up is schematically shown in 
Fig.JIJa). The SQUIDs and gate are controlled by a com- 
bination of dc and rf signals, mixed together by bias tees. 
The rf signals are generated by synchronized arbitrary 
waveform generators (AWG in the Figure) , guided to the 
sample by coaxial lines, attenuated and thermalized at 
different temperature stages. Low-pass filters (LP) with 
60 dB attenuation up to 40 GHz are placed at the sample 
stage. The dc wiring consists of 160 f2 surface-mount re- 
sistors, 2m-long lossy coaxial lines, and lm-long twisted 
pairs. At room temperature, a voltage bias VJ, produced 
by a floating digital-to-analog converter (DAC) is applied 
through a divider, and current is read out by a transcon- 
ductance amplifier with sensitivity 10~ 10 A/V. The sam- 
ple is protected by two nested rf-tight shields in order to 
prevent microwave irradiation from higher-temperature 
stages. 

We measure the sluice in the supercurrent branch, 



close to Vb — 0. We apply control pulses at frequency 
/ (/ = 80 MHz unless otherwise stated). V is nulled 
by minimizing the current flowing without applying any 
pulses, thus compensating for the voltage induced by the 
current amplifier. Application of the flux pulses alone 
was found not to shift the zero-bias point, nor induce 
any additional current. This rules out the presence of 
rectification effects, which could have been introduced, 
e. g., by capacitive couplings between the on-chip coils 
and the leads. To further reduce the influence of voltage 
fluctuations, the pumped current I p is detected as the 
difference between the currents measured while pumping 
in opposite directions. 

Evidence of single Cooper-pair pumping is presented 
in Fig.|2j where the pumped charge per cycle Q p — l v j f 
is plotted versus the offset n° g and the peak-to-peak am- 
plitude An g of the gate drive. Since the pumping cycle 
of Fig. [2] is fully adiabatic, we expect Q p to be quan- 
tized in units of Cooper pairs, the first correction being 
of the order of J mini J max [11] . From energy-diagram 
considerations, it is easy to show that the regions of 
constant Q p are diamond-shaped in the offset-amplitude 
plane. The regions are delimited by the family of curves 
A n s — 2 1 rig + m\, where m can be any integer. 

A remarkable feature of Fig.[2]is that Q p is 2e-periodic 
in the gate charge (that is, the size of the diamonds is 1 
unit along the n g and 2 units along the An g axis). Pre- 
vious measurements with the Cooper-pair sluice [T^l [T7] 
reported le-periodic plateaus in the pumped charge plot- 
ted against the amplitude of the gate drive (no depen- 
dence on the gate offset was reported). The authors as- 
cribed the observed periodicity to quasiparticle poison- 
ing. Quasiparticle poisoning [18] has been intensively 
studied in systems closely related to the sluice, the single 
Cooper-pair transistor [T5H2T] and the single Cooper-pair 
box [221 [23] ■ All these devices feature a superconduct- 
ing island in the Coulomb-blockade regime. For a given 
position of the gate, there are two metastable states for 
the island ("odd state" and "even state"), differing by the 
presence of one quasiparticle. Noncquilibrium quasiparti- 
cles generated in the leads (at temperatures T -C A/ks, 
the thermal population of quasiparticle states is unim- 
portant) drive transitions between the two states ("par- 
ity fluctuations"), shifting the gate charge by exactly half 
a Cooper pair. In the sluice, this results in le-periodic 
pumped charge, provided that the time scale of parity 
fluctuations is intermediate between the pumping period 
and the acquisition time [13 EZ]. In our case, the clean 2e 
periodicity observed implies that the device is not "poi- 
soned" by quasiparticles. We ascribe this improvement 
to two factors: efficient microwave shielding and clean 
grounding of the probe leads. The importance of mi- 
crowave shielding has been emphasized in other recent 
experimental works [231 [23] . 

While the data of Fig. [2] are clearly not affected by 
nonequilibrium quasiparticles, the latter may come into 
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FIG. 3. (Color online) Crossover between pure Cooper-pair 
and mixed Cooper-pair-quasiparticle dynamics, (a) Pumped 
charge Q p versus gate offset n° g for increasing gate amplitudes 
An g (bottom to top), (b) Circles: same data as in (a), pro- 
jected on the rig axis and plotted versus An g . Dashed line: 
asymptotic adiabatic-limit expectation. In both panels, the 
first three pumping plateaus are indicated by arrows. 



play as a result of nonadiabatic pumping. One way of 
making the pumping nonadiabatic is to increase the am- 
plitude of the gate modulation. This increases both the 
effective speed of the drive (continuously) , and the num- 
ber of tunnelling events involved (discretely). We do it 
in Fig. [3|a), where Q p is plotted versus w? g with An g 
taking a series of values in the range of 0.1 and 7. The 
same data are also plotted against An g in Fig. [3Tb). The 



data show a clear crossover between pure Cooper-pair 
and mixed Cooper-pair-quasiparticle dynamics. Up to 
about An g — 3, Q p is 2e-periodic in n®, as in Fig. |2] 
The pumping plateaus are also 2e-periodic in An g . They 
show up as straight lines in Fig. [3ja), and nodes in 
Fig. |3jb). The first three plateaus are indicated by ar- 
rows. The crossover takes place between about An g = 3 
and An g = 5, where the pattern is blurred. Finally, for 
An g > 5 a clear periodicity is restored, but the period has 
doubled. These data show that quasiparticle poisoning, 
while initially absent, can be induced by noncquilibrium 
quasiparticles generated by a nonadiabatic drive. The 
link between loss of adiabaticity and quasiparticle poison- 
ing is strengthened by the fact that the crossover accom- 
panied by a reduction in Q p with respect to the adiabatic- 
limit expectation, as shown in Fig. [3^b) for An g > 5. 
As nonadiabatic transitions occur in the sluice as missed 
Cooper-pair tunnellings, they leave a detectable trace in 
overall magnitude of the pumped charge. 

We understand the generation of noncquilibrium quasi- 
particles as a multi-step process: LZT to an excited 
state followed by relaxation via Cooper-pair breaking and 
quasiparticle tunnelling. When the drive is nonadiabatic, 
one or more Cooper-pairs may fail to tunnel as dictated 
by the gate, leaving the island in an excited state. Even 
if there are no quasiparticles in the leads near the junc- 
tions, quasiparticle tunnelling from/into the island is still 
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possible provided the energy gain AiS* of the process [55] 
exceeds the energy cost 2A required to break a Cooper 
pair. In our case, due to the moderate charging energy 
of our device (Eq — 0.33A) this is possible when the 
charge occupation of the island differs from that of the 
ground state by at least two Cooper-pairs. After the tun- 
nelling of a quasiparticle, the island is left in a metastable 
state with an odd quasiparticle number, which may then 
decay by tunnelling of a second quasiparticle. Overall, 
this mechanism is similar to the well-known Josephson- 
quasiparticle cycle [57] , with the nonadiabatic gate drive 
playing the role of an effective voltage bias in dynamically 
creating nonequilibrium. 

We now take a closer look at individual LZTs. To do 
so, we choose to pump a single Cooper pair at a time. 
This reduces the dynamics to that of a two-level system. 
Transitions between the adiabatic ground (g) and excited 
state (e) may occur due to Landau-Zcner tunnelling at 
avoided crossings [see Fig. Qa)]. For a single crossing 
the transition probability Plz = e 
the adiabatic parameter S = J^ ax /Hv 
Josephson coupling of the active SQUID at the cross- 
ing and v = Ecdrig/dt the rate of change of the energy 
difference between diabatic (charge) states. The full tun- 
ability of our device gives us the possiblity to control the 
degree of adiabaticity in several independent ways. In 
we plot Q p versus J max , for the case n g = 
\n g — 0.45. The traces are taken at different pump- 
ing frequencies in the range of 70 and 120 MHz, normal- 
ized to the asymptotic pumped charge Qq, and vertically 
offset for clarity. All traces approach a constant value for 

large values of J max , and monotonically decrease to for 
/ / jo v n 

'Jmax/ •-'max ~ u - 

In order to understand the data quantitatively, 
we must consider the interplay between LZTs and 
environment-induced decoherence. Decoherence effects 
in the sluice have been studied theoretically with a 
master-equation approach [THl |55J HH] ■ Numerical calcu- 
lations indicate that even in the presence of modest deco- 
herence, the driven Cooper-pair sluice approaches a qua- 
sistationary state [55] after a few tens of cycles. Since our 
acquisition time extends over about 10 6 pumping cycles, 
we identify Q p with the quasistationary charge pumped 
by the sluice in the presence of its electromagnetic envi- 
ronment, averaged over a great number of cycles. 

Let us now discuss the role of decoherence in deter- 
mining the quasistationary state and hence Q p . In the 
present set-up, the sluice is directly connected to room- 
temperature leads. Hence we expect a fast dephasing 
time, of the order of a few ns. This implies that sub- 
sequent LZTs are totally uncorrelated, which rules out 
the possibility of Landau-Zener-Stiickelberg interference 
|30] . The latter might instead play a role at higher pump- 
ing frequencies and/or in a closed- loop geometry [14] , 
Besides dephasing, we also have to consider relaxation, 
which can take place via environment-assisted (inelas- 
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FIG. 4. (Color online) Nonadiabatic pumping and Landau- 
Zener transitions, (a) Instantaneous energies of the ground 
and first excited state versus time for a pumping period. 
In the model described in the text, nonadiabatic transitions 
are localized at level crossings, and occur with probability 
Plz- Decoherence induces complete dephasing between sub- 
sequent transitions, and relaxation via inelastic Cooper-pair 
tunnelling (both are indicated by wavy arrows), (b) Normal- 
ized pumped charge Q p /Qo versus Josephson coupling J max , 
for a set of frequencies in the range of 70 and 120 MHz. The 
experimental traces (squares) are vertically offset by 0.2, and 
plotted together with the best-fit of Eq.[I]to them (solid lines). 
(c,d) Asymptotic pumped charge Qo versus / (c) and A pa- 
rameter versus 1/f (d) for three different measurement sets 
(squares, circles, and triangles). The solid line in (d) is a fit 
to the data of the expression X = fo/ /• 



tic) Cooper-pair tunnelling through the SQUIDs. The 
energy scales and base temperature for our system are 
such that the transition rates to the excited state are 
exponentially suppressed, so that the environment is ef- 
fectively at zero temperature. As a result, relaxation 
tends to keep the sluice in the ground state [TB], coun- 
teracting LZTs. In general, determining the explicit de- 
pendence of Q p on Plz for an arbitrary environment is 
a formidable task, which can only be undertaken using 
numerical methods. In order to provide a simple analyt- 
ical formula to be compared to our data, we shall fur- 
ther assume that (i) relaxation fully takes place between 
subsequent anti-crossings, so that the system approaches 



every anti-crossing in the ground state, and (ii) it takes 
place at the two anti-crossings equally, so that on average 
it does not contribute to the pumped charge. Under (i) 
and (ii), the problem greatly simplifies, and we find that 
Q p /2e = 1-P LZ . 

Based on this model, in Fig. |4ja) we fit to each exper- 
imental trace the expression 



Q p (x) = Qo [l - exp (-27rAa; 2 )] 



(1) 



where x — J max I Jmax ano - A is a free parameter. The 
excellent agreement between the data and our single- 
parameter fit provides strong evidence that the departure 
from the adiabatic limit takes place via LZTs, and that 
our understanding of decoherence effects, albeit simpli- 
fied, is essentially correct. 

The dependence of Qo and A on / is shown in 
Fig. |4jc,d). Different symbols refer to three different 
measurement sets, taken on the same sample using dif- 
ferent generators and/or during different cooldowns. Al- 
together, they span the frequency range of 40 to 120 
MHz. This range is limited from below by acceptable 
signal-to-noise ratio, and from above by the sampling 
rate of our rf generators. As the three sets yield con- 
sistent results, we discuss them together. The behavior 
of Qo [Fig. |4^c)] is not completely clear; it may stem 
from frequency-dependent attenuation in the lines. By 
contrast, A [Fig. |3Jd)] displays a clear trend, a steady 
increase with decreasing /. The solid line in Fig. |4jd) 
is a best-fit to the data of the expression A = fo/f, 
yielding / — 180 MHz. Our model predicts fg — 
{J max ) 2 /(AhaE c /±.n g ), where a ~ 4.3 is a parameter 
proportional to the steepness of the gate rise. Using our 
estimated Eq and J^ ax , we obtain 120 MHz. We 

find the agreement between f and /q satisfactory. The 
discrepancy between the two corresponds to an effective 
Josephson energy about 20% smaller than its nominal 
value. Among other factors, this may well be due to in- 
fluence of noise/decoherence not captured by our simple 
model. 

Our results demonstrate the feasiblity of using Cooper- 
pair pumps to investigate coherent effects in driven quan- 
tum systems, and pave the way for further observations. 
By embedding the sluice in a fully superconducting loop 
(as e.g. in Refs. Q2] and and increasing the pump- 
ing frequency up to a few hundred MHz, we expect to 
achieve coherence times extending over several pumping 
periods. This would allow challeging proposals such as 
Landau-Zener interferometry with geometric phases [14] . 
characterization of decoherence induced by an engineered 
environment |28j . and measurement of the Lamb shift 
[32], to be readily implemented. 
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